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We derive here the equation of state for quark matter with a nontrivial vac- 
uum structure in QCD at finite temperature and baryon density. Using thermofield 
dynamics, the parameters of thermal vacuum and the gluon condensate function 
are determined through minimisation of the thermodynamic potential, along with 
a self-consistent determination of the effective gluon and quark masses. The scale 
parameter for the gluon condensates is related to the SVZ parameter in the context 
of QCD sum rules at zero temperature. With inclusion of quarks in the thermal 
vacuum the critical temperature at which the gluon condensate vanishes decreases as 
compared to that containing only gluons. At zero temperature, we similarly obtain 
the critical baryon density for the same to be about 0.36 fm -3 . 
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I. INTRODUCTION 



Hot dense hadronic matter in relativistic heavy ion collisions or in early universe is 
likely to behave as quark gluon plasma. The study of this in the context of quantum 
chromodynamics (QCD) is nontrivial and nonperturbative [|l]||. For gauge fields only we had 
seen earlier || that such a vacuum structure with gluon condensates can emerge dynamically. 
We examine here the same problem including quarks at zero or finite temperature, and at 
finite baryon density, but again with only gluon condensates . We expect that the effect 
of gluons will dominate vacuum structure due to the colour factor. 

The method used here is variational and nonperturbative, since only equal time algebra 
is taken as input. The ansatz functions for variation are determined with the minimisation 
of energy density at zero temperature, and, of the thermodynamic potential at finite tem- 
peratures and baryon densities. This is a continuation of the earlier programme || using 
thermofield dynamics method (TFD) Q to study the problem at finite temperature and 
leads to a simple generalisation of Bogoliubov transformations as a part of the technology 
at zero temperature. 

We organise the paper as follows. In section II, we consider quark matter at finite 
temperature using as stated the method of thermofield dynamics [f|]. In section III, we 
calculate the thermodynamic potential for the system of quark matter at finite temperature 
and baryon density and then minimise the same to derive the equation of state. In section 
IV we discuss the results. Similar techniques with explicit ground state construction have 
also been applied to study hot nuclear matter || where scalar isoscalar two pion condensates 
replace the u-field and for chiral symmetry breaking || with quark condensates. 

Nonperturbative nature of QCD has been studied for a long time in the context of 
condensates |1],@J7]] or scale symmetry breaking || as well as effective field theories |6|,[| or 



lattice gauge theories 10|. The present approach is complementary to the above. 
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II. QUARK MATTER AT FINITE TEMPERATURE 



We shall very briefly recapitulate the notations of Ref. [Q], now with quark fields so that 
finite baryon density effects can be included. Let us consider the QCD Lagrangian given as 



where 



£ £ gauge ^matter ^inti 



(1) 



c 



gauge 



l -G a ^(d^W\ - d v W% + gf ahc W\W c u ) + \g\ v G^\ 



(2a) 



matter 



(2b) 



and 



An* = gH l Y w ft> 



(2c) 



where W a M are the SU(3) colour gauge fields. We shall quantise in Coulomb gauge and 
write the electric field G a oi in terms of the transverse and longitudinal parts as 



r~ia 1 r~ia , o ra 

<J Oi — ^01 + °il j 



where f a is to be determined. We take at time t=0 
W a i{x) = (2tt)- 3 / 2 

and 



dk -> -> t 

(a a i(k) +a a i(-k) )exp(ik.x) 

2u(k) 



\| ^(-« a *(£) + a a i(-k) j ) exp(ik.x) 



T G a 0l (x) = (2ir)- 3 / 2 t I dk 
where, u(k) is arbitrary ]ll| and for equal time algebra we have 



a^y^Y] = 5 ab Aij(k)5(k - k'), 



(3) 



(4a) 



(4b) 



(5) 



with 
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Ay® = % - ^- 
The equal time quantization condition for the fermionic sector is given as 



We now also have the field expansion for fermion field if) at time t=0 given as 

1 



(6) 



(7) 



ip\x) 



(2vr) 3 / 2 

where U and V are given by [|12 



U r (k)c\ r (k) + V s (-k)5 l Is (-k) e zk - x dk, 



(8) 



U r (k) 



cos 



X(k) 



\ 



a ■ k sin j 



Ulr] V s (-k) 



a.ksm ^p- 



\ 



cos 



x(fc) 



Via 



(9) 



/ 



where the function x{k) could be arbitrary. Here we approximate the same as cosx(^) = 
rriQ/e(k) and sinx(A;) = k/e(k), with e(k) = (k 2 + nig) 1 / 2 . For free fields e.g., itiq = tuq. 
However, for interacting fields, we shall determine iriq in a self consistent manner as will be 
discussed later. The above are consistent with the equal time anticommutation conditions 
provided fl2 ] 



Ir 



(£),4(#) t l+ = s rs s ij 5(k-k') = [4(^),4(^) t ]- 



(10) 



where i and j refer to the colour and flavour indices . 

In Coulomb gauge, the expression for the Hamiltonian density, T 00 from equation (1) is 
given as 
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T 



00 



T G a 0i T G a 0t + -W a t (-X7 2 )W a i + gf abc W?Wpm 



+ ^f abc r f w?w<w t e wf + ^rmn 

+ ipi-ifdi + m )ip - gilyf—W^ :, 



(11) 



where : : denotes the normal ordering with respect to the perturbative vacuum, say | vac >, 
defined through a a i{k) \ vac >= 0, c l Ir {k) \ vac >= and c l Ir (kY \ vac >= 0. In order to 
solve for the operator f a , we first note that 



f 1 = ~W\ - g f abc (V 2 f\w\ d^o). (12) 

Proceeding as earlier [RJ with a mean field type of approximation we obtain, 

V 2 W a (x)+g 2 f abc f cde < vac', (3 | WM^d^y^W^x) I vac',(3> d 3 W\{x)) 

= J a o(x), (13) 

where, 

JS = Qf^wfGl, - gjrf (14) 

We note that at zero temperature, | vac';f3 = oo >=| vac' > was the nonpertubative 
ground state as discussed in Ref. 0. The extra fermionic contributions in eqautions (11) 
and (14) may be noted, and thus the expressions of Ref. M will get modified at finite 
temperatures and densities. We define | vac' > through a unitary transformation, in a 



similar manner to Gross-Neveu model considered earlier |T3[, given as 



| vac >= U | vac >, (15) 

where 

U = exp(fi f - B), (16) 

In Ref. 0, it was shown that at zero temperature, we may have 

B ] = \j f($a\$a\(-k)*dk, (17) 

where f(k) describes gluon condensates. For the consideration of vacuum destabilisation we 
should also take quark condensates |J. However, we shall not do the same here because of 
the following reason. Our method here shall consist of minimisation of the thermodynamic 
potential. We had earlier determined the temperature dependent mass like term for the 
gluon in a self consistent manner 0. We shall here also need to determine the same for 
the quark fields. The numerical computation with a variation for the vacuum structure of 



gluon and quark condensates becomes forbidding. Besides this, we believe that the gluons 
being in the adjoint representation shall contribute more strongly to QCD interactions, and 
therefore dominate the vacuum structure. However, we do include the effect of quark sector 
for vacuum structure at finite temperatures through thermofield dynamics [|J. 

In the following, we shall consider the effect of temperature as well as finite baryon density 
on the behaviour of the gluon condensates for vacuum structure. We use the method of 
thermofield dynamics to consider the above which is convenient for our purpose while dealing 
with operators and expectation values. The thermal average of an operator is replaced by 
expectation value in an extended Hilbert space associated with thermal doubling [[| . For the 
present case of including temperature and finite baryon density effects, the thermal vacuum 
is given as 

\vac',(3 >= U G (/3)U Q (/3)\vac' > (18) 

where U G and Uq are unitary operators involving thermal excitations of gluons and quarks 
respectively. For the gluon sector, we have the old expression g 

U G {/3)=e X p(B G {/3? -B G {/3)), (19) 

with 

B G (f3) ] = J 9{k,f3)b a i (kfb a i {-kfdk. (20) 
In addition, for quark sector we have, 

U Q {/3)=e X p{B Q {/3? -B Q {/3)), (21) 

with 

b q {(3) ] = J [e-dpyu^dj^-ky + 9 + (k,p)^ Ir (kyc i Ir (-k)^ dt (22) 

In the above the underlined operators || correspond to the extra Hilbert space in TFD. 
Further, 9(k,/3), 9±(k,f3) are arbitrary functions to be determined from the minimisation of 
the thermodynamic potential. For example, for free fields these functions are given by 
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sinh 2 9(k,P) 



exp{/3u(k, /3)) - 1 



(23a) 



and 



sin 2 9±{k,P) 



1 



exp(/3(e(A;,/3)± //)) + ! 



(23b) 



where cu(k,P) = \J k 2 + tuq, e(k, (3) = ^Jk 2 + m 2 Q and fi is the quark chemical potential. 
However, for interacting fields these will be different and we shall approximately determine 
them in a self consistent manner. We have seen that such a structure as in equation flT8| ) in- 
troduces a thermal Bogoliubov transformation for gluon fields 0. We now have, in addition, 
the parallel transformation in quark sector given as 



V 



cVr(fc) 

c l Ir (—k) 
c l Ir (k) 



\ 



( 



cos sin Q_ 
- sin cos 6- 



\ 
















cos 9 + sin # H 

— sin 9+ cos 9. 



\ 



( 



c'lriKP) 

c l Ir (—k, P) 

c} r (-k,P) 
~d Ir (k,pf 



(24) 



Our job now is to evaluate the expectation value of T 00 with respect to | vac; P >. For 
this purpose, we have the earlier eqautions [[| 

< vac'; p \: W a i {x)W b j {y) :| vac'; P >= 5 ab x (2tt)- 3 f dke^^ ^y^y Ay(jfe), (25) 



< vac';P\ : T G^(x) T G b 



r 0j yyj . | Vac' ; P > 

5 ab x (2tt)- 3 / dke ik - ( - £ -®A ij (k)u(k,P)F_(k,P). 



(26) 



In the above the temperature dependant F±(k,P) are given as 

F ± (k,P) = cosh2fl(sinh 2 /(fc) ± Smh ^ fc) ) + sinh 2 9{k,p) 



(27) 



where sinh 2 9(k, P) is given by equation (|23a|) . For the quark fields we have the parallel 
eqautions given as 



< >vac>,P= (27T)- 3 ^" 



A_(k,P)) ^-^"^dk, 



(28a) 



< 



where 



A±(k, 0) = =f^ (sin 2 6L - sin 2 # + ) + (7°cosx + a.k sin x)( sin 2 6L + sin 2 6> 4 

Using equations (|TT|), fl2"E| ) and (28), we then obtain the expectation value of T 0( 
respect to | vac; (3 > as 

e (/3) = < vac'; |: T 00 :| vac'; /3 > 

= C F {(3) + C x (/3) + C 2 {(3) + C 3 (/3) 2 + CM, 

where 

Cj?(/3) = <: 4>(-iYdi + m )ip :> mc ',/3 

6 f k^ elk — * _» 

y (sin 2 M fc > 0) + sin2 M fc > (3))(m Q m + k 2 ), 



and 



tt 2 7 e(Jfe) 

I 

2 

4t / uj(k)k 2 F_(k,f3) dk 
% 2 J 



Ci(f3) - <: - T G a 0i T G a oi :> vac >.p 
4 



C 2 ((3) = <: I^^-V 2 )^, :> mc ^ 



c 3 ((3f = <: ^ 2 r fec r e/ ^iy c ,4y e 4 ^ :> vac , p 

(2g f k 



,2 \2 



ir 2 J u(k,p) 



F + {k,f3) dk 
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In the above, 



G 1 (k,P) = 3g 2 / dqF + (\q\,/3) F_(\k + q\, (3) 



w{\k + q\,(3) 

w(|g|,/9) 



X 



(1+ 



{q 2 + k.g} 2 
q 2 (k + q) 2 



(32a) 



G 2 (k,(3) = -g 2 Jdq (l + 



+ 



q.(q- k) 



e{q,(3)e{q-k,(3) e(q, (3)t{q - k, (3) 
x ( sin 2 (q, (3) sin 2 Q_ (q -k,/3) + sin 2 6 + (q, (3) sin 2 6 + (q - k, (3) 

m Q ((3) 2 q.(q- k) 



- 1 



e(g, /?)e(g - fc, /?) e(g, /3)e(g - fc, /3) 
(sin 2 0-(£/3) sin 2 0+(<f- £,/?) + sin 2 9+(q,p) sin 2 fc,/3)) 



(32b) 



and 



<f>(k, (3) 



dk ' F +{ k,p)(k 2 + k--^:S )2 ^ k+k ' 



2kk' 



k-k' 



(32c) 



]tt 2 J co(k',P) 

The expressions above are the same || as earlier, except for the additional contribution of 
Cp{(3) from the fermionic terms, as well as the fermionic contribution through G2{k,f3) in 
the expansion for C±{f3) arising from the auxiliary fields. 

Since for interacting fields, the form of the functions u(k,/3) and e(k,(3) are not known, 
we parametrise them in the free field form with temperature dependent effective mass pa- 
rameters for the gluon and quark fields given as 



oo(k,{3) = yJk* + m G (P)>; e(k,P) = yJk* + m Q (/3¥ 



(33) 



with ma(f3) and mg(/3) are to be calculated self consistently as below. 

We identify the gluon mass m G (j3) as earlier P,|T4|| through the self consistency 
requirement that 

2g 2 f k 2 



u(k,/3) 

This is derived through single contraction contribution from the quartic terms 



(34) 
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We similarly identify the effective quark mass, mg from the sum of the single contractions 
of the term quartic in the field ip of T 00 given by equation (|11|) as well as from the mass 
term mo'0'0- Writing the sum of the single contraction terms of T 00 i nt as 

C(f)M a/3 (v,)^(^), (35) 

we identify the effective quark mass, as m,Q{(3) = m +mg(/3), with uiq given through 

the relation 

MapiVx) , - , n = m' Q {P)6 a p. (36) 
Writing ijj l {x) in the momentum space as 

V>*(x) = (2tt)- 3/2 / ft(k)exp(ik.x)dk } (37) 



we thus have 

4 # 2 /" dkdk'dkx 



^j:\^°A + (k^))jm 



3 (2tt) 6 7 ( fel _ A; /)2 + 0( A . 1 _ A ./ j/ g) 
(2tt)- 3 / ^ a {k)M aP {k')^{k')dkdk' . (38) 



Hence, with the identification M a/3 (k') ^ ^ = mg(/3)5 Qi g, the effective quark mass, mq(f3) 
as given by 

m«(« = »„ - x ^g(^ M £,/») + -n»M*»- (39) 

This is the se// consistency requirement for the quark mass, and is solved iteratively where 
the input rriQ(P) of the right hand side through e(fc, /?) becomes equal to the output mQ(/3) 
of the left hand side. 

We shall extremise over the thermodynamic potential containing eo(/3). For this purpose, 
as earlier HI we shall take 



sin 



hf(k) = Ae~ Bk2/ \ (40) 



which corresponds to taking a gaussian distribution for perturbative gluons in nonperturba- 
tive vacuum ||. The energy density, e (/3) in terms of the dimensionless quantities x = yfBk, 
jj,G = \fBm G ((3) , jiQ = yBmQ(P) and y = -4= then gets parametrised as 
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e (A, (3) = W2 {I F {y) + h{A, y) + I 2 (A, y) + I 3 (A, yf + h(A, y)) 



where 



If(v) 



6 r x 2 dx 



7T 



e(x,y) 



>in 2 0_ (x, y) + sin 2 6 + (x, y)){n Q {y)no + x 2 



and 



/ < (A,,) = 4x(2 T )-/jg Gl( Vit Gj( f ,) . 

j x^ + (p{x,y) 



In the above, Gi(x,y) = G(x,y) of 0, and 

Ml/) 2 



G 2 (x,y) = -g / c/x 



1 + 



+ 



x'.(x' — x) 



X 



e(f, y)e(f - f , y) e(f, y)e(£' - x , y) 
( sin 2 9-(x , y) sin 2 ^-(af — x, + sin 2 6+(x, y) sin 2 0+(a? — x,y 



- 1 



e(f , y)e(f - f, y) e(x' , y)e(x' - f, y) 
( sin 2 6_(x', y) sin 2 6* + (x' — x, y) + sin 2 9 + (x', y) sin 2 6L(af — x, y) 



(41) 



(42a) 



(42b) 



(43) 



The other expressions I\(A,y), I 2 (A, y), I 3 (A,y) and <fi(x,y) are the same as in M. The 
above integrals contain Hciy) an d Hq(3j) which are determined from the self consistency 
requirements |§ 



VG(y) 2 = % I 4^ ( A " e ~ x2 + Ae^(l + A 2 e-*y/ 2 ) ( 1 + 
ti z J u [x, y V /V 



w(a;,2/) 
1 



exp(ya;(x,j/)) - 1, 



exp(yo;(x,2/)) - 1 



(44) 



and 



r 



dx- 



x ^^-( sin 2 64x,y) + sin 2 6 + (x,y)). (45) 



37r 2 J x 2 + (f)(x } y) e(x, y) 
Hg{v) and Hq{i/) are solved through an iterative procedure B. We shall now calculate the 



thermodynamic potential and minimise the same. 
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III. EXTREMIS ATION OF THERMODYNAMIC POTENTIAL AND RESULTS 



At zero temperature, we had considered extremisation of energy density to obtain the 
vacuum structure. At finite temperatures and baryon densities, the relevant quantity for 
extremisation is the thermodynamic potential, which at temprature T = 1//3 is given as 0] 

T{A, (3) = e (A, (3) -^{S G + S F ) - hbPb- (46) 

Here ps = 3/x is the baryon chemical poential corresponding to the baryon number density 
Pb given as, with two quark flavours, 

1 1 ( 

PB = -x2x3x2x — — J ( sin 2 6_ - sin 2 6 + )dk, (47) 



and p as in equation (|23b| ) is the quark chemical potential. Further, in the above, €q(A, (3) 



is as given in equation (|30|) , and the entropy densities S G and Sp for the gluon and quark 
fields are given as 

S G = -2 x 8 x (2tt)- 3 J dk[ sinh 2 6 log(smh 2 6) - cosh 2 6 log(cosh 2 6)) (48) 

and 



S F = -3x2x2x ( 27r ) J dk[sin 2 e.(k, (3)log(sm 2 6_(k, (3)) 

+ (cos 2 6- (k, (3))log(cos 2 6^ (k, 13)) + sin 2 9+(k, f3)log( sin 2 9+(k, /?)) 

+ (cos 2 6 + (k,f3))log(cos 2 6 + (k,[3))). (49) 



Clearly, the factor 2 x 8 in fl48 ) above comes from the transverse and colour degrees of 
freedom for the gluon fields and the factor 3 x 2 x 2 in (|4^) comes from the colour, flavour 
and spin degrees of freedom for the quarks as well as for the antiquarks. For the minimisation 
of thermodymic potential we may scale out the dimensional parameter l/I? 2 and write 

1 



T{A^) = — F 1 (A,y) 



E 1 
1 

52 



F(A, y) - -(S G (A, y) + S F (A, y)) - p' B p' B 

y 
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(50) 



where F(A,y) = B 2 eo(A, (3), p' B = y/BjiB, p' B = B 3 ^ 2 p B and the entropy densities Sc(A,y) 
and Sf(A, y) in dimensionless units are given as 

1 , 1 



off 1 

Sa(A,y) = — J / x 2 dx\ ( — — -)log( 

ix l J L exp \yuj\x, y)) — 1 



and 



exp (yu(x,y))-l exp (yu (x, y ) ) - 1 
^ + exp (yuj(x,y)) - l^ l ° 9 ^ + eX p (yu(x,y)) - 1^} ^"^ 



tt 2 / X ^{^exp (y(e(x,y) - p) + l^° 9 ^exp (y(e(x,y) - fj) + 1^ 
+ (! / / 7 \ rTT)^^ " ' 



exp (y(e(x, y)-[i) + l exp (y(e(x, y)-fi) + l' 

+ ( — , , , \ . \ , J M- ' 



exp (y(e(x, y) + /i) + 1 exp (y(e(a;, y) + /i) + 1 

+ 0- m \m ui )/o ^ (1 FT7 ~TT n ■ j j- ( 52 ) 

exp (y(e(x, y) + n) + 1 exp [y{e{x, y) + fi) + l J 

We note that for each A, the gluon and quark masses are determined self consistently 
through equations (44) and (45) for the evaluation of the right hand side of equation (50) . 
We extremise Fi (A, y) of equation fl5T?D with respect to the parameter A and obtain the 
optimum value of A as A min at a given temperature T and baryon density p#. In fig.l we 
plot A min versus T at zero baryon density for y 2 /47r = .8 against the earlier curve which 
was without quark excitations in the thermal vacuum. We find that the critical temperature 
is about 205 MeV in contrast to 275 MeV of Ref. where quark excitations were not 
included. We have taken here the coupling constant y 2 /47r = 0.8 and the Lagrangian 
quark mass m =300 MeV as a typical value for the constituent quark mass. Such a result is 
similar to that of lattice QCD where inclusion of two flavours of quarks decreases the critical 
temperature [15| to about 150 MeV [16j from about 235 MeV of quenched approximation 



|T7|| . For quark mass m less than 300 MeV, Tq decreases below 205 MeV, but for very 
small quark masses, the numerical calculations tend to be unreliable. A min is plotted in 
fig.2 as a function of baryon density ps for temperatures and 100 MeV. We find that A min 
decreases with increase in baryon density p#, and vanishes at and above a critical value, 
(ps)crit of pb- The values of (ps) cr it for temperatures and 100 MeV are 0.36 fm~ 3 and 
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0.88 fm~ 3 respectively. Clearly in fig. 2 A m i n decreases with ps and hence in fig. 1 for finite 
Pb the A min versus T curve will always lie below the solid curve, yielding a smaller critical 
temperature. In fig. 3, we plot the effective gluon mass as a function of pb for the above 
temperatures. We see that it starts decreasing with increase in the baryon density, ps, and 
for zero temperature, it becomes zero at and above (pB) C rit=0.36 fm -3 and for T=100 MeV, 
it decreases and remains almost a constant above (ps) C rit=0.88 fm~ 3 . This is because, when 
the condensate function vanishes, the effective gluon mass depends only upon the thermal 
contributions as in equation (44) and is independent of the baryon density. We then plot in 
fig.4, the effective quark mass, rriQ as function of pB, which starts decreasing with increase 
in p B and there is a discontinuity at {psjcrit above which the quark mass decreases slowly. 
This discontinuity in the mass is associated with abrupt vanishing of condensate function 
with A min approaching zero at the critical baryon density. Similar behaviour was observed 
for the effective nucleon mass by Ellis et al in Ref. ||, where however the critical baryon 
density at which the gluon condensates vanishes was much higher, being of the order of 
2.37 fm~ 3 . This discontinuity might indicate a possible link between vanishing of gluon 
condensates and chiral symmetry restoration. 

We next estimate the SVZ parameter using the value of A min . This at any temperature 
is given as 



with I\(A, y), I^A, y), Iz{A, y) and Ii(A, y) as in equation (41) ||. In fig. 5, we plot the SVZ 
parameter given above as a function of ps- It starts decreasing with increase in ps and for 
T=0, becomes minimum at {pb)ctU and then increases with increase in ps- It may be noted 
that although the condensate function vanishes at (pB)crit, the SVZ parameter is nonzero 
and increases with density. This is because the contribution from the fermionic sector to 
I4 is negative and increases in magnitude with density. Similar behaviour is also seen for 
T=100 MeV except that the magnitude here is smaller as we have extra positive thermal 
contributions to I\ and I4. 




vac 



B 2 7T 2 L 




h(A, y) + I 2 (A, y) + I 3 (A, yf - I 4 {A, y) 



J A=A, 



(53) 
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We also calculate the pressure P given as [ITS 



P(P) = -F(A,P) 



A. — A-n 



(54) 



In fig. 6, we plot the pressure as function of the baryon number density, pg for the tem- 
peratures and 100 MeV. We see that the pressure increases with increase in pg. We may 
compare the present nonperturbative results with the perturbative estimation of the pressure 
|i~9|l given as 



Pert 2?r2 

_ ^1 

7T' 



Pfkf(p 2 f — 2.5m 2 ) + 1.5miln 



1.5 x pjkf — mJn( 



m f 



k 4 



(55) 



where k 2 = p 2 — m 2 and pf = /is/3 is the baryon chemical potential. In the above we 
have taken a s =0.8 and considered the case for two flavours of quarks with masses 300 MeV 
each. This perturbative equation of state for zero temperature is plotted as the dashed curve 
in Fig. 6. As may be seen, the present equation of state is stiffer than the perturbative 
equation of state. 



IV. DISCUSSIONS 

In the present paper, we have extended the study of nontrivial ground state structure 
of QCD with gluon condensates at zero and finite temperatures to the case of nonzero 
baryon densities. The calculations as earlier are done in Coulomb gauge. We would have 
also liked to include the quark condensates for the vacuum structures in these calculations. 
It would have been nice to demonstrate that this effect is not large, as we believe to be the 
case. However, simultaneous consideration of gluon condensates and quark condensates with 
appropriate self consistency requirements becomes computationally prohibitive and has not 
been attempted here. 

The modified gluon condensate function and the quark distribution functions are ob- 
tained here through minimisation of the thermodynamic potential with, as stated, a self 
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consistent determination for the quark and gluon mass parameters. At different tempera- 
tures, the gluon condensates disappear at critical values of the baryon density. Inclusion 
of thermal excitations in quark sector as here decreases the critical temperature from 275 
MeV to 205 MeV even for zero baryon density, and shall be lower at finite baryon densities. 
At zero temperature, the above critical baryon density is about 0.36/fm 3 . The results are 
similar to that of lattice QCD [^5l ~ |l7| . 

We note that the present approach is based on QCD Lagrangian of quarks and gluons 
without the introduction of effective fields |§, but with a nonperturbative variational ap- 
proach and an explicit vacuum structure 0. This may be relevant for quark gluon plasma 
where baryon structures are likely to dissolve at high densities or temperatures |20"|-p2]j . 
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Figure Captions 

Fig.l: We plot A min as function of T in MeV . The dashed curve corresponds to the case 
of not including thermal excitations in the quark sector. Lowering of critical temperature 
with inclusion of quarks in the thermal vacuum may be noted. 

Fig. 2: We plot A min as function of ps in fm -3 for temperatures and 100 MeV. 

Fig. 3: We plot the gluon mass, me in MeV as function of ps in fm -3 for temperatures 
and 100 MeV. 

Fig. 4: We plot the fermion mass, mq in MeV as function of ps in fm -3 for temperatures 
and 100 MeV. The discontinuity in the same at critical baryon density may be noted. 

Fig. 5: We plot the SVZ parameter in units of 10~ 2 GeV A as function of ps in fm -3 for 
temperatures and 100 MeV. 

Fig. 6: We plot the pressure, P in units of MeV/fm 3 as function of pb in fm™ 3 for tem- 
peratures and 100 MeV. The perturbative equation of state as given by equation ( |55D for 
temperature, T=0 is given as the dashed curve of the same figure. 
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